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±^ Abstract 

g 

I By means of theoretical modeling and experimental synthesis and characterization, 

?H we investigate the structural properties of amorphous Zr-Si-C Two chemical compo- 

, '■^ , sitions are selected, Zro.31Sio.29Co.40 and Zro.60Sio.33Co. 07- The amorphous structures 

1-H are generated in the theoretical part of our work, by the stochastic quenching (SQ) 

> 

\0 method, and detailed comparison is made as regards structure and density of the 

OC experimentally synthesized films. These films are analyzed experimentally using X- 

cr^ ray absorption spectroscopy, transmission electron microscopy and X-ray diffraction. 

O 

CN Our results demonstrate for the first time a remarkable agreement between theory 

^ and experiment concerning bond distances and atomic coordination of this complex 



C^ 



/\^ amorphous metal carbide. The demonstrated power of the SQ method opens up 

avenues for theoretical predictions of amorphous materials in general. 
Keywords: Amorphous materials, Ab initio calculations. X-ray absorption 
spectroscopy. Metal carbide glasses. Atomic structure 
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1. Introduction 

Deposition of thin films from the vapor phase may lead to amorphous structures. 
A group of materials which often exhibit non-crystalline structure is various tran- 
sition metal-metalloid films such as Cr-C, Zr-Si, and Fe-B.[l]-|3] Interestingly, they 
are also important components in metallic glasses formed during rapid quenching 
from a melt, suggesting that the tendency to form amorphous structure is related 
to inherent properties of the elements and interactions between these elements. To 
obtain a more detailed understanding of the formation of amorphous structures and 
glasses we need to develop theoretical methods to model the structure and bonding 
in these non-crystalline environments. This can be achieved by e.g. molecular dy- 
namics (MD) |3], Monte Carlo simulation [5], and reverse Monte Carlo simulation. [6] 
Unfortunately, these methods suffer from being computationally expensive or to rely 
on interatomic potentials with sometimes questionable accuracy. In this study we 
investigate another approach, namely the stochastic quenching (SQ) method, [3 |8| 
which combines computational efficiency and accuracy in describing the chemical 
interaction. In this method the atoms are placed randomly in the calculation cell, 
and then they are relaxed using first-principles density functional (DFT) calcula- 
tions until the force on every atom is negligible. It has been shown that this way of 
generating amorphous structures is possible. [3 |H] However, the present investigation 
is the first in which a detailed comparison between experiment and the SQ theory 
is made in terms of structural properties, such as bond length and nearest neighbor 
coordination, in particular for a material with complex chemical interactions. 

In this study, we have selected the Zr-Si-C system as an amorphous model sys- 
tem, since transition metal carbides and silicides are known to have complex chem- 
istry, with competing metallic and covalent bonds. |9] We deposit amorphous films 
using magnetron sputtering from elemental targets and confirm the structure exper- 
imentally with several techniques. The choice of elements is based on several facts. 
First, Zr is a well-known base element in many metallic glasses. The atom has a 
large radius (1.59 A), favorable for amorphous structure formation when combined 
with several other elements with smaller radii. [TD] Second, the Zr-Si system is well- 
known to produce amorphous thin films with magnetron sputtering for potential 
use as e.g. diffusion barriers. [3j The addition of carbon to a ternary Zr-Si-C film 



will create a multicomponent system with different atomic radii favorable for glass 
formation. [10] Also, in a ternary Zr-Si-C film a wide variety of bond types are formed 
(metallic Zr-Zr, and covalent Zr-C, Zr-Si and Si-C bonds), which makes it possible 
to create a network structure which should further favor an amorphous structure. 
The potential for applications in Zr-Si-C has recently been discussed. [TT] 

The paper is organized as follows. In Section |2j we describe the theoretical 
methods used in this study, and summarize the most important details of the calcu- 
lations. The experimental methods are described in Section [3] The theoretical and 
experimental results are presented in Section |4} 

2. Theory 

The amorphous structures were generated in the theory by the stochastic quench- 
ing (SQ) method. This method is designed to describe amorphous structures in gen- 
eral, although there is less experience with it for complex materials with competing 
natures of the chemical binding. The amorphous structures are obtained by means 
of the following two steps of the SQ method: Firstly, the atoms are placed randomly 
in the calculation cell under the constraint that no pair of atoms are closer than a 
small value (typically 0.4 A). This constraint is required in order to avoid numerical 
problems in the first few steps in the calculation. Secondly, the positions are relaxed 
by means of a conjugate gradient method until the force on every atom is negligible. 
We found that the self averaging of 150-200 atoms usually properly describes the 
electronic properties at the thermodynamic limit. |H1 fT2HT6] In the present study 
we used 200 atoms for most of the calculations, and made a few calculations with 
400 atoms, for comparison. 

The first-principles calculations were performed by means of the projector aug- 
mented wave [171 [18] method as implemented in the Vienna ab initio simulation 
package (VASP). [ISHH]. This method is based on DFT. [221 123] The exchange- 
correlation energy was calculated using the generalized gradient approximation with 
the Perdew, Burke, and Ernzerhof functional. [24j The calculations were considered 
converged when the potential energy difference between atomic iterations was less 
than 10~^ eV/atom, and the forces on each atom were typically less than 0.005 
eV/A. A plane-wave energy cutoff of 400 eV was employed. The calculations were 
performed using only the F k-point. 



3. Experimental 

The Zr-Si-C film was deposited by non-reactive, unbalanced, DC magnetron 
sputtering in an ultra-high vacuum chamber (base pressure of 1*10~^ Pa). Separate 
2 inch elemental targets (Kurt J. Lesker Ltd) with a purity of 99.999% for Si and 
C and 99.2-99.7% (grade 702) for Zr were used. The Ar plasma was generated at a 
constant pressure of 0.4 Pa with a deposition temperature of approximately 350 °C 
and bias of -50 V. All depositions were made using a rotating substrate holder to 
ensure a homogenous composition. Prior to deposition the Si and Si02 substrates 
were cleaned using ultrasonic bath in 2-propanol and ethanol for 5 minutes each 
and dried with nitrogen gas. The deposition rate was 40-70 A/min depending on 
composition and the films for XRD, XPS and resistivity analysis was deposited 
with a thickness of ~0.3 yum. The chemical composition was analyzed by XPS 
using a Physical Systems Quantum 2000 spectrometer with monochromatic Al Ka 
radiation and using sensitivity factors determined from binary reference samples of 
known composition. The crystallinity of the film was studied using XRD and TEM. 
The XRD analysis was made using a Siemens D5000 diffractometer and Cu Ka 
radiation. A gracing incidence (GI) scan with an incidence angle of 1° was used to 
increase the signal from the thin film in relation to the substrate. For top view TEM 
the sample was prepared by depositing 40 nm thin film on a NaCl substrate. The 
NaCl was then dissolved in distilled water after deposition releasing the deposited 
film, which then was put on a carbon-coated copper grid and analyzed using a JEOL 
JEM-2100 operated at 200 kV. For the resistivity measurements each sample was 
measured six times and mean value and 95% confidence interval was calculated. 

The XAS measurements were performed at BM26A (Dutch-Belgian beamline 
"DUBBLE") |25] of the European Synchrotron Radiation Facility (ESRF) in Greno- 
ble, France. The incident energy was selected using the <111> refiection from a 
double Si crystal monochromator with bendable 2nd crystal for sagital (horizontal) 
focusing. Vertical focusing and rejection of higher harmonics was achieved with Si 
mirror with Pt layer at an angle of 1.8 mrad relative to the incident beam. The 
incident X-ray beam had a fiux of approximately 1 x 10^^ photons/s on the sample 
position. The XAS data were measured in fiuorescence mode by 9-element mono- 
lithic Ge detector [26j at Zr K edge (17988 eV) at room temperature. The Zr-Si-C 



thin films on Si02 substrates were conventionally mounted on a thin tape. The en- 
ergy calibration was performed on Zr foil. The data were recorded with a constant 
0.05 A^^ step in fc-space, in the energy range 17.75-18.85 keV. The intensity was 
normalised to the incident flux. The total energy resolution was estimated to be 
-3.5 eV. 

Spectra were analysed using the iXAFS implementation of IFEFFIT, Athena, 
and Artemis. f27\ fIE\ Model data were obtained from averages of many calculations 
of EXAFS data obtained from local structures present in the description of the 
amorphous material, using the FEFF software. |29]43T] 

4. Results and Discussion 

On the experimental side, amorphous structures are typically generated by fast 
cooling from melts or with quenching of atoms during growth directly from the gas 
phase. Two films with composition Zro.60Sio.33Co.07 and Zro.31Sio.29Co.40 (according 
to our XPS analysis) were here deposited using magnetron sputtering. In our XRD 
measurements both films exhibited identical types of structure. Figure [l] shows 
the GI-XRD diffractograms of the Zro.60Sio.33Co.07 and Zro.31Sio.29Co. 40 films. Two 
broad features, one at 36° and the other one at 60° is observed in the diffractogram. 
Another small peak is observed at 61° for the Zro.60Sio.33Co. 07 film, which is only a 
contribution from the Si(lll) substrate. This type of diffractograms with only broad 
features is characteristic for an X-ray amorphous material, where these features are a 
consequence of the close-range order in the material. Fig. |2]shows the TEM image of 
the Zro.60Sio.33Co.07 film. No indication of crystallinity is visible and the film shows a 
completely amorphous structure. Also, the energy filtered electron diffraction (inset) 
shows only broad, featureless rings without indication of any crystalline contribution. 
TEM and XRD together confirm that both films are amorphous. 

The amorphous structures were generated in the theory by the SQ method. We 
calculated the energies of 50 stochastic structures of both compositions studied in the 
experiments, at different volumes, with supercells of 200 atoms (Fig. |3]). We found 
an amorphous structure in all generated configurations. No partial crystallization or 
porosity was observed (Fig. |4]). The calculated energies at constant V show roughly 
a Gaussian distribution for both compositions (insets in Fig. |3]). We have determined 
the equilibrium volume (Vq) and bulk modulus (Bq) by fitting an exponential Morse- 



type function to the average energies of Fig. pi and obtained Vo=13.72 A^/atom 
for Zro.31Sio.29Co. 40, and Vo=19.08 A^/atom for Zro.60Sio.33Co.07- Accordingly, our 
calculated equilibrium densities are f)o=4.99 g/cm^ for Zr0.31Si0.29C0.40, and go=5.64: 
g/cm^ for Zr0.60Si0.33C0. 07- This is in line with the experimental densities of the corre- 
sponding binary crystalline phases: ZrC (6.730 g/cm'^ L32J), ZrSi2 (4.883 g/cm^ [33]). 
and ZrsSis (5.998 g/cm'^ I31])- We calculate -Bo=86 GPa for Zro.31Sio.29Co.40, and 
i?o=116 GPa for Zro.60Sio.33Co.07- Within accuracy, Bq in Zro.60Sio.33Co. 07 is about 
the same as that in crystalline ZrSi2 (-BzrSi2=114.8 GPa [35]), Bq in Zro.31Sio.29Co.40, 
however, is noticeably lower than that. Bq of both amorphous materials is smaller 
than that of crystalline Zr5Si3 (i?zr5Si3 = 152 GPa [34J). As expected, the theoretical 
B of amorphous Zr-Si-C is significantly smaller than that of the crystalline carbides 
(Szrc=223 GPa, 5sic=211 GPa). 

We note that in amorphous materials, in contrast to crystalline materials, not 
only bond lengths can change upon compression, but also bond angles may vary, 
even to the extent that atoms may rearrange under compression. This invariably 
reduces the curvature of the energy versus volume relationship, when compared to a 
crystalline material with a reduced bulk modulus as the result. It is possible that this 
explains the much lower bulk modulus for amorphous Zro.31Sio.29Co.40 (-Bo=86 GPa) 
when compared to ZrC. The presence of Si also changes the bonding situation which 
may contribute to the relatively low bulk modulus of Zro.31Sio.29Co.40- It should also 
be mentioned here that the larger bulk modulus of Zro.60Sio.33Co.07 {Bq=116 GPa) 
compared to that of Zro.31Sio.29Co.40 is somewhat surprising. The compound with 
more carbon has a higher density, and is expected to have more of strong Si-C bonds, 
which is expected to increase the bulk modulus, compared to a material with lower 
C content. Our calculations give the opposite trend. As we show in our study, the 
structures generated by the SQ method agree very well with observations. Hence 
it remains to be seen if the calculated trend in bulk modulus is a true effect, or 
possibly the result of the numerical treatment of the 50 configurations considered 
at each volume. 

In order to compare the structural properties of the theoretical and experimental 
amorphous phases, in Fig. [5] and [6] we compare the simulated X-ray absorption fine 
structure (EXAFS at the Zr K-edge) signals of Zro.31Sio.29Co.40 and Zro.60Sio.33Co. 07 
to the experimental results. The simulations of the EXAFS spectra were performed 



by means of real space multiple-scattering theory as implemented in the FEFF code 
|36j . The size of our cells prohibited the use of self-consistent potentials for the 
multiple scattering calculation of the X-ray absorption coefficients so we distributed 
overlapping atomic potentials on the sites in our amorphous 200 and 400 atom cells. 
The multiple scattering path expansion was used as the full multiple scattering 
approach loses accuracy at high k-values. In order to simulate periodic boundary 
conditions, we constructed a 3x3x3 supercell and calculated the normalized EXAFS 
signal x?(k) from each Zr atom in the center cell. The index i denotes individual 
Zr atoms, and the index C denotes individual 200 atom or 400 atom cells obtained 
by means of the SQ technique. 

For each cell C, the calculated normalized EXAFS spectra from each individual 
Zr atom x?(k) were averaged to form the final x*"(k)- The theoretical spectra were 
calculated at K and at room temperature by means of the correlated Debye model. 
The Debye temperature for the amorphous material is unknown, so we estimated a 
value of 600 K for the calculations that was close to the Debye temperatures of ZrC 
(614 K - 680 K [371138]), ZrSis (495 K [39]), and ZrgSig (480 K \M\). 

The theoretical curves of the 50 stochastic structures with 200 atoms (grey lines 
in Fig. |5]and[6]) are very close to each other, especially in the small k region, and with 
increasing k they become more diffuse. We find that the overall agreement between 
the averages of the 50 individual structures (red lines) and the experimental signals 
(dashed black line) is good. In particular, we obtain excellent agreement between 
theory and experiment in the small k region. When comparing Fig. |5]and Fig. |6]we 
note that Fig. |6] reproduces observations almost perfectly, whereas a disagreement 
may be seen in Fig. |5J in the region of 6.5-7 A^^. We made a few additional 
calculations with 400 atoms per unit cell. As expected, due to the better statistics, 
the curves of the individual stochastic structures are closer to each other than for the 
200 atom calculations [Fig. [s] . The average of the individual structures practically 
does not change compared to the 200 atom results, which clearly shows that a 200 
atom unit cell is sufficient to describe amorphous Zr-Si-C 

Returning to the disagreement between theory and experiment in the region of 
6.5-7 A~^, for the sample with low Zr concentration, it can be analyzed further by 
making a Fourier transformation of the data in Fig. Isl This transform (data not 
shown) shows that two bond lengths at distance just above 2.5 A are located at 0.2 



A shorter distances in the theoretical data compared to the experiment. Although 
this may be viewed as a smaller error to an overall rather satisfactory agreement 
between experiment and theory, it can most likely be traced back to the very fast 
quenching of atomic coordinates in the SQ method, which is a somewhat slower 
process in the experimental samples. This may be realized by noting that for the 
low Zr concentration case, the two peaks for which the SQ method makes a 10 % 
error are associated to the average Zr-Zr and Zr-Si bond lengths. In this sample 
we expect few nearest neighbor Zr-Zr coordination, and a reduced amount of Zr-Si 
bonds, due to a preferred Zr-C bonding in combination with a rather large amount 
of C in the sample. However, an inspection of the simulated geometries gives that 
the SQ method does find a finite amount of Zr-Zr nearest neighbor coordination, due 
to the rapid quenching. For the sample with larger Zr concentration the amount of 
C is very low, and hence neither simulations nor experiment can exclude a nearest 
neighbor Zr-Zr coordination, resulting in a better agreement between theory and 
experiment for the EXAFS data. 

To demonstrate that the agreement between theory and experiment is really due 
to an amorphous structure, we also calculated the EXAFS spectra of some known 
crystalline compounds of Zr, Si, and C (Fig. [T]). In this EXAFS calculation we 
used the experimental crystal structures. The Debye temperatures used for these 
crystals were 680 K (ZrC) [3H], 495 K (ZrSig) [SOj, 480 K (ZrsSig) \M\, and 600 K 
(Random alloy) j^. The simulated EXAFS curves of these compunds do not agree 
at all with experiment. We also examined a case, where the calculated composition 
(Zro.4oSio.2oCo.4o) differs slightly from the experimental one (Zro.43Sio.30Co.27) to a 
small extent. For such a composition, the simulated curve does not agree with the 
experimental one, neither in the small, nor in the large k region (data not shown). 
This shows that both accurate structural information as well as precise chemical 
composition is needed to reproduce the experimental EXAFS data. The data in 
Figs. |5] and |6] are the first direct comparison on a structural level, between the SQ 
method and observations, and the agreement shows that a computationally efficient 
theory provides accurate atomic coordinates for a complex amorphous metal carbide. 

In the following, we present additional structural details about amorphous Zr- 
Si-C The theoretical average partial radial distribution functions (RDF) calculated 
for amorphous Zro.31Sio.29Co.40 and Zro.60Sio.33Co.07 (Figs, [sjandjo]) show short range 



order up to 6 A for the C bonds (Zr-C, Si-C, and C-C), and up to 7-8 A for the longer 
Zr-Zr, Zr-Si, and Si-Si bonds. In Table [l] the calculated average nearest neighbor 
distances are compared to distances in typical crystalline phases. Atoms are con- 
sidered nearest neighbors (and also counted in the coordination numbers later), if 
they are closer than the first minimun in the corresponding partial RDF curves. 
We used the following cutoff distances for Zro.31Sio.29Co.40 (and Zro.60Sio.33Co.07): 
4,zr-zr=4.11 (4.29) A, 4,zr-si=3.60 (3.61) A, (/,,zr-c=3.00 (3.15) A, 4,si-Si=3.52 
(3.10) A, 4,si-c=2.33 (2.41) A, and 4,c-c=l-89 (1.92) A. For most pairs the dis- 
tances are comparable to the corresponding crystalline bonds. The pairs with shorter 
bonds correlate with a reduced coordination of carbon. A lowered coordination 
means that the remaining bonds are stronger than in the crystalline phase. Such 
bond-shortening mechanism was observed in MAX phases. [H] 

On the average, Zr atoms are coordinated by 15.4 atoms in Zro.31Sio.29Co.40, and 
by 14.9 atoms in Zr0.60Si0.33C0.07 (Table |2|. Both are larger than the coordination 
number in hexagonal close-packed (hep) structure (?7,hcp=12). Such high coordina- 
tion is not unusual in amorphous structures. |l2] Zr atoms are mostly coordinated by 



Zr atoms (Tableland Fig. 10a), which may be primarily due to their high concen- 



tration. Si atoms have a slight preference to have Zr neighbours (Fig. 10b), which 
is pronounced in Zro.60Sio.33Co.07, where the Zr-Si coordination is more than double 
that of the Si-Si coordination (Table |2|, and Zr has almost ten times larger number of 
Si nearest neighbors than Si (see Table [3]). Due to the small C content of amorphous 
Zro.60Sio.33Co.07, nzT,c, ""-Si.c and nc,c (defined in caption of Table [2]) are well below 



one, and C atoms are mostly surrounded by Zr atoms (Fig. 10c). As expected, the 
number of C-C, Si-C and Zr-C bonds is significantly larger in Zro.31Sio.29Co. 40 (Table 
[3|. Here C atoms are mostly coordinated by Zr atoms, and on the average, each of 



them has at least one C nearest neighbors. Figure 10 shows that the structure can 
be described as a network structure of metal-metal, metal-carbide and metal-silicon 
bonds. 

To examine the electronic structure of amorphous Zr-Si-C, we calculated their 



electronic density of states (DOS). Figure 11 shows the calculated the electronic 
density of states (DOS) of the representative structure of Zro.60Sio.33Co. 07 at the 
equilibrium volume. The lowest lying states at around -12 eV originate from C 2s 
states, and mainly Si 3s orbitals contribute to states between -10.5 and -6 eV. The 



large band between -5 and 0.6 eV is built up by Zr Ad, Si 3p, and C 2p states. 

We calculate a large DOS at the Fermi level [Ep) for both compositions, i.e. that 
amorphous Zr-Si-C has a metallic character. This is in line with experimental ob- 
servations: we measured 3.20±5 /if2m resistivity for Zro.60Sio.33Co. 07, and 10.30±18 
/ifim for Zro.31Sio.29Co.40- The states at E-p are dominated by Zr Ad orbitals, with 
some contribution from Si 3p and C 2p orbitals, i.e. amorphous Zr-Si-C is conduct- 
ing, and the Zr d band has the most important role in electrical transport. In a 
simple approximation, the density of states at the Fermi level, N{Ep), is directly re- 
lated to the electrical conductivity. We calculate N{Ep)=0.785 states/eV/atom for 
amorphous Zro.60Sio.33Co.075 and A^(-Ef)=0.353 states/eV/atom for Zro.31Sio.29Co.40, 
suggesting that Zro.60Sio.33Co. 07 has a higher conductivity. Indeed, we measured 
lower resistivity for Zro.60Sio.33Co. 07 than for Zro.31Sio.29Co.40- We obtain higher DOS 
at Ep in amorphous Zro.60Sio.33Co. 07 than that in crystalline MAX phases, e.g. in 
TiaSiC (0.36 states/eV/atom [H]), TigGeC (0.43states/eV/atom [43j), ZrsInC (0.30 
states/eV/atom |ii]). 

5. Conclusions 

By means of the SQ method, we have generated structures for amorphous 
Zr0.31Si0.29C0.40 and Zro.60Sio.33Co.07, and measured the EXAFS spectra, TEM and 
XRD diffractograms of samples with the same chemical composition. TEM and 
XRD confirm that both compositions are amorphous. We show that the calculated 
EXAFS spectrum is in excellent agreement with experiments, which serves as proof 
that theory reproduces accurately the structural properties of amorphous materials. 
This demonstrates for the first time that the SQ method provides reliable atomic 
coordinates of amorphous materials, even for very complex ternary systems, with 
competing chemical interactions, involving metallic and covalent bond-formation. [9] 
Our findings open up a new avenue to perform fast and accurate theoretical simu- 
lations for a wide class of amorphous materials. 
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Figure captions 

Figure 1: (Color online.) X-ray diffraction patterns for Zr-Si-C films with 7 at.% C 
(lower) and 40 at.% carbon (upper). 

Figure 2: (Color online.) TEM micrographs of Zr-Si-C film with 7 at.% carbon. 

Figure 3: (Color online.) Theoretical average energies (-Eav) calculated for 
Zro.31Sio.29Co.40 (upper) and Zro.60Sio.33Co.07 (lower), as a function of the volume, 
V. In the insets the energy distribution of 50 stochastic configurations are shown 
at Vo, where the energies were smeared by normalized Gaussian distributions with 
standard deviation of 0.01 eV. 

Figure 4: (Color online.) One stochastic structure of amorphous Zro.60Sio.33Co.07, 
calculated at the equilibrium volume, having a total energy closest to the average 
energy of 50 structures of the same volume (representative stucture). Zr, Si, and C 
atoms are displayed by silver, yellow, and cyan, respectively. For clear representation 
of the structure, bonds are not displayed in the figure. 

Figure 5: (Color online.) The simulated (grey lines) and experimental (dashed lines) 
EXAFS signals of amorphous Zro.31Sio.29Co.40 (in arbitrary units). The averages over 
50 individual configurations for the 200 atom case (upper) and over 5 configurations 
for the 400 atom case (lower) are shown by red lines. 

Figure 6: (Color online.) The simulated (grey lines) and experimental (dashed line) 
EXAFS signals of amorphous Zro.60Sio.33Co. 07 (in arbitrary units). The average of 
all simulated cells are shown by the red line. 

Figure 7: (Color online.) Simulated EXAFS signals (in arbitrary units) for some 
crystalline compounds compared to the measured spectrum of amorphous 
Zro.31Sio.29Co.40- The crystalline alloy is a random substitutional alloy with the same 
composition as the amorphous alloy distributed on a NaCl-type crystal lattice with 
lattice constant a=4.698 A. 

Figure 8: (Color online.) Average partial radial distribution functions calculated for 
amorphous Zro.31Sio.29Co.40- The supercell size in the calculation was 200 atoms. 

Figure 9: (Color online.) Average partial radial distribution functions calculated for 
amorphous Zro.60Sio.33Co. 07- The supercell size in the calculation was 200 atoms. 



Figure 10: (Color online.) Examples of calculated local structures around Zr (a), 
Si (b), and C (c) atoms in amorphous Zro.60Sio.33Co.07- Zr, Si, and C atoms are 
displayed by silver, yellow, and cyan, respectively. 

Figure 11: (Color online.) Electronic density of states (in arbitrary units), calculated 
for the representative structure of amorphous Zro.60Sio.33Co. 07 at the equilibrium 
volume. The top panel shows the total DOS of the 200 atom structure. The partial 
s, p, and d DOS curves are averages calculated for the different atom types. The 
Fermi level (Ep) is denoted by vertical dashed line. 



Table 1: Average nearest neighbor distances (c^Typei-Typeii in A) calculated for amor- 
phous Zro.31Sio.29Co.40 and Zro.60Sio.33Co. 07- For comparison, the distances in the 
corresponding crystalline materials are also listed. 





Zro.31Sio.29Co.40 


Zro.60Sio.33Co.07 


Crystalline 


Zr-Zr 




3.22 




3.18 


3.179° 


Zr-Si 




2.75 




2.75 


2.724^ 2.723^= 


Zr-C 




2.32 




2.28 


2.349^ 


Si-Si 




2.41 




2.45 


2.352^ 


Si-C 




1.86 




1.89 


1.887^ 


C-C 




1.41 




1.41 


1.545^ 1.42'* 



°: in Zr of hexagonal P63/mmc structure [32] 

'': in ZrSi2 of Cmcm structure |l5] 

'^■. in ZrsSis of PGs/mcm structure 

'^: in ZrC of NaCl structure [32] 

^: in Si of diamond structure [32] 

^■. in SiC of F43m structure [32] 

^: in diamond [32] 

^■. in graphene [32] 



Table 2: Calculated average coordination numbers of 50 stochastic structures for 
amorphous Zro.31Sio.29Co. 40 and Zro.60Sio.33Co. 07- Type! atoms (in rows) are coordi- 
nated by nxypei.Typeii uumber of Typell atoms (in columns) on the average. 





Zro.siSio.^ 


J9C0.40 


Zro.ec 


Sio.33Co.07 




Zr 


Si 


C 


Zr 


Si C 


Zr 


6.6 


4.6 


4.2 


10.2 


4.1 0.6 


Si 


4.9 


3.6 


1.6 


7.4 


1.5 0.2 


C 


3.3 


1.1 


1.1 


5.2 


0.9 0.1 



Table 3: Calculated minimal (TTimm)) maximal (mmax' 



and average (m 



average , 



number 



of nearest neighbors in amorphous Zro.31Sio.29Co.40 and Zro.60Sio.33Co.07- Numbers in 
parentheses give the percentages of the total number of nearest neighbors, where 
the minimal and maximal values are calculated for the corresponding individual 
stochastic structures. 





TTijniii 


Zr3iSi29C 

^^max 


40 
^^averagc 


TTljiiin 


Zr6oSi33C 

^^max 


7 

^^ average 


Zr-Zr 


191 
(19.4) 


217 

(22.2) 


204.4 
(20.7) 


597 
(48.5) 


623 

(50.7) 


612.7 

(49.4) 


Zr-Si 


258 
(26.1) 


305 

(30.8) 


285.1 
(28.8) 


467 
(38.0) 


511 
(41.0) 


490.2 
(39.6) 


Zr-C 


235 

(23.9) 


282 
(27.9) 


260.7 

(26.4) 


59 

(4.8) 


78 
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72.0 

(5.8) 


Si-Si 
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120 

(12.3) 
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(10.5) 


38 
(3.1) 


64 
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50.2 
(4.1) 


Si-C 


77 
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108 
(10.8) 


90.7 
(9.2) 
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(0.6) 


20 
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13.2 

(1.1) 


C-C 
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(5.9) 


44.0 
(4.5) 
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